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Abstract Measurements of the ionized Ca n K line are one of the major 
resources for long-term studies of solar and stellar activity. They also play a 
critical role in many studies related to solar irradiance variability, particularly 
as a ground-based proxy to model the solar ultraviolet flux variation that may 
influence the Earth's climate. Full disk images of the Sun in Ca II K have been 
available from various observatories for more than 100 years and latter synoptic 
Sun-as-a-star observations in Ca n K began in the early 1970s. One of these 
instruments, the Integrated Sunlight Spectrometer (ISS) has been in operation 
at Kitt Peak (Arizona) since late 2006. The ISS takes daily observations of solar 
spectra in nine spectra bands, including the Ca n K and H lines. We describe 
recent improvements in data reduction of Ca n K observations, and present time 
variations of nine parameters derived from the profile of this spectral line. 
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1. Introduction 

The diagnostic value of Ca II K measurements for the study of solar magnetism 
and irradiance variability is well established (e.g., IFoukal et al.l 120091 and ref- 
erences therein). These measurements provide information on energy input to 
solar chromosphere and its long-term variation. The intensity of the K-linc is a 
proxy for total magnetic flux fe.g-. IOrtiz and RastU2005]) . and it can be used to 
represent solar variability in (extreme) ultraviolet, (E)UV (e.g.. :Foukal | 12002)) . 
Solar observations are also used to interpret and calibrate stellar measurements 
of Ca II K and H spectra in terms of stellar activity and cycles fe.g.. lRuttenlll9841 



Giampapa et al. 120061 Vieytes, Mauas, and Dia z 2009; IBuccino et al.i 120111 
Grunhut et q/.ll2T)10p . As shown in some studies (e.g. IDonahue and Keill IT995; 
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|Hasler, Riidiger, and Staude| I2002[) , the connection between the modulation of 
the disk-integrated Ca n K -line flux and solar surface differential rotation may 
open the possibility of using these measurements to determine the differential 
rotation of other main-sequence stars. 

In a broader context, Ca II K observations have been used to estimate dis- 
tances to early-type stars ( Megier et al.\ I2009P , and as diagnostics of velocity 
fields in red giant stars ( |Vieytes et al. l2~01ip . It has also been suggested to use 
recurring chromosphcric emission visible in Ca II K as a tool to detect massive 
extra-solar planets orbiting very close to their central stars QLenz, Reiners, and KiirsterU 
I20TU1) . 

Since 1974, a long-term program of monitoring the sun-integrated Ca n K line 
has been operating at the National Solar Observatory at Kitt Peak ( |Livingston, Milkey, an d Slaug hter] 
H977l|Livingston et q/.|l200"7]) and since 1977 at Sacramento Peak (jKeil and Wordenl | 
IT9"g4"l |Keil, Henry, and Fleck] UM^i . Starting December 2006, disk-integrated 
sunlight spectra in Ca n K line are recorded on a daily basis by the Integrated 
Sunlight Spectrometer (ISS). The ISS is one of three instruments, with the Vector 
Spectro-Magnetograph (VSM) and the Full Disk Patrol (FDP), comprising the 
Synoptic Optical Long-term Investigations of the Sun (SOLIS) - a synoptic facil- 
ity for solar observations operating at NSO/Kitt Peak ( |Keller, Harvey, and Giampapa| [ 
12003a| |Keller, Harvey, and The Solis Team| [2003Ell. A Solar-Stellar Spectrograph! 
project started at Lowell Observatory also provides regular observation of the 
solar K line since 1994 (IHall and Lockwoodl [T9981 |Hall, Lockwood, and Skiff] 
l2T)07l) . 

Currently the ISS instrument provides measurements in nine different spectral 
bands, including the Ca II K band centered at 393.37 nm. |Livingston et al\ 
(201(5]) have shown that the ISS Ca K index compares well with both the Ca 
II K McMath-Pierce spectral line scans and the 10.7 cm radio flux. Here we 
provide the first detailed description of the data reduction and describe recent 
improvements in the determination of the nine different parameter time series 
derived from the Ca II K spectra. These improvements result from changing 
most of the original ISS code that processes the observed Ca n K line profiles 
and then computes the parameters. The improvement is more significant for 
some of the parameters, while others benefit only marginally from these changes. 
Sections 2 and 3 provide a brief description of the ISS instrument, the observing 
sequence, and the data reduction. In Section 4, we describe the derivation of 
nine parameters representing properties of Ca II K line profiles, and compare 
the results of improved data analysis with previous data reduction. Section 5 
provides a brief summary of the results of this study. 



2. The ISS Instrument 

The Integrated Sunlight Spectrometer(ISS) operating at the National Solar Ob- 
servatory at Kitt Peak (Arizona) is designed to obtain high spectral resolution (R 
= 300,000) observations of the Sun as a star in a broad range of wavelengths of 
350 nm -1100 nm. Beginning December 1, 2006, the ISS takes daily observations 
in nine spectral bands centered at CN head band 388.40 nm, Ca n H (396.85 
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Table 1. Spectral bands measured by the ISS instrument. Each spectral band is centered 
at Ao, with a bandwidth given by AA. The average linear dispersion over time, dX/dx, is 
also given. The start date (fifth column) indicates the first day ISS began to observe in 
that particular spectral band. 



Type Aq, nm AA, nm d\/dx, pm/pixel Start date 



CN band 


388.40 


0.58 


0.564 


December 4, 2006 


Ca II K 


393.37 


0.55 


0.541 


December 1, 2006 


Ca II H 


396.85 


0.53 


0.522 


December 4, 2006 


C I 


538.00 


0.84 


0.824 


December 4, 2006 


C I (with iodine lines) 


538.00 


0.84 


0.824 


January 7, 2008 


Mn I 


539.41 


0.83 


0.816 


December 4, 2006 


H-alpha 


656.30 


1.14 


1.12 


August 31, 2007 


Ca II 


854.19 


1.61 


1.58 


December 13, 2006 


He I 


1083.02 


1.65 


1.61 


December 4, 2006 


NaDl 


589.59 


0.98 


0.956 


March 23, 2011 



nm), Ca n K (393.37 nm), C I 538.00, Mn I 539.41 nm, H-alpha 656.30 nm, 
Ca II 854.19 nm, He I 1083.02 nm, and NaDl (589.59 nm). In addition, high 
resolution spectra formed by the sunlight passing through an iodine vapor cell 
were recorded beginning January 7, 2008. Iodine spectra are over imposed on the 
solar spectral range centered at C I 538.00 nm. Table [T] lists currently observed 
spectral lines. 

Observations of integrated sunlight with the ISS are accomplished through the 
use of a fiber optic feed. A small optical system (a lens of 8-mm diameter) in- 
stalled on the side of main mount of SOLIS/VSM focuses a 400 micron diameter 
image of the Sun on the input face of a 600 fi diameter fiber. The fiber assembly 
transmits light to a McPhcrson 2-m Czerny- Turner double-pass spectrograph 
located in a temperature-controlled room below the telescope. The output beam 
from the fiber consists of sunlight scrambled in both angle and position so 
that any output angle from any position on the exit face is well-integrated. 
A prism pre-disperser isolates the desired wavelength band. The spectrograph 
employs a 316 g/mm grating blazed at 63.5 degrees and a movable, back- 
illuminated 512 x 1024 CCD as the focal plane detector. Additional information, 
including transmission characteristics of the fiber optics, the instrument's block 
diagram, mechanical and optical layout are available from the SOLIS web site at 
http:/ /solis. nso.edu/ISSOverview. html Daily observations can be accessed from 



|http://solis. nso.edu/iss 



3. Data Processing 

For each spectral band, observations are taken in four different positions by 
moving the CCD in respect to the spectra in the spectrograph focal plane: the 
first image (image 1, 512 x 1024 pixels in size) is taken with the CCD centered 
at a corresponding wavelength (for a given spectral band), image 2 is shifted by 
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1 pixel in the direction of dispersion, image 3 is shifted by 11 pixels, and image 4 
by 129 pixels. It takes about 5 minutes to complete the full cycle of observations 
for most spectral bands; for Ca n K observations it takes about 12 minutes, 
and for He I 1083.0 nm about 15 minutes. Four images are used both for the 
purpose of flat fielding and to create the final spectrum for a given spectral band. 
Because in four images the CCD is shifted in respect to fixed spectral features, 
the same pixels are exposed to a different level of light. Assuming that the solar 
source remains constant during the observing cycle (3 - 15 minutes), the images 
can be used to derive flatfield and CCD gain based on an algorithm described in 
|Toussaint, Harvey, and Toussaint|(|2003p . A dark field is taken for each observed 
spectral range at the beginning of four-image cycle by closing the intermediate 
shutter of the double pass system. The same dark field is used for all four images 
taken in each spectral band. 

After applying dark and flat fields, spectral features on four spectral images 
are aligned, and the astigmatic images are averaged in the spatial direction 
to create the final spectral line profile for a given spectral band. In this aver- 
aging, each spatial spectrum of the four 512 x 1024 pixel images is weighted 
by the square root of its total intensity. The final products are intensity as 
a function of pixel number (spectral line profile) and its error for each pixel 
(computed as standard deviation) . At the next step of data reduction, the pixel 
position of selected spectral lines is determined. Pixel positions are converted 
to nanometers using known wavelengths of selected spectral lines from high 
resolution spectra taken with the NSO Fourier Transform Spectrometer (FTS, 
| Wallace, Hinkle, and Livingston) I2007p . 

Spectra reduced to this level still show a slight gradient in intensity in wave- 
length direction. For spectral bands that include continuum, this gradient is 
removed by a linear fit to the continuum and scaling the spectra by the fit. For 
the Ca II K and H lines, ISS observations do not include continuum. These two 
spectral bands are scaled using intensities at two narrow spectral bands situated 
in blue and red wings of these lines. Mean intensities in these two bands are 
scaled to match spectral line profiles taken by the FTS. For Ca II K spectra 
these two bands are: 393.147-393.153 nm in the blue wing, and 393.480-393.500 
nm in the red wing. The intensity profile is then scaled so that the mean residual 
intensity between 393.4958-393.5134 nm is 0.1756. 

Reduced spectra are saved in FITS format as double array with three axes: 
wavelength in Angstroms, continuum normalized intensity, and error in intensity 
for each pixel. These level 1 FITS files constitute the basic ISS data product. 



4. The ISS Ca n K Parameter Time Series 



The importance of Ca n K parameter time series for monitoring solar and stellar 
activity has been established by many decades of studies (e.g. Wilson and Vainu Bappufl 
119571 IWilsoni 119681 White and L ivingston I1978p . Currently, nine parameters are 
routinely derived from each measured 393.37 nm ISS spectrum. The quality of 
these parameter-time series depends strongly on our ability to measure the linear 
dispersion of the instrument over time. 
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Table 2. Observed line profiles within the ISS Ca K band and cor- 
responding wavelengths (A) used in the determination of the linear 
dispersion. 



Element 


A, nm 


Element 


A, nm 


Element A, nm 


Fe I 


393.11251 


Fe I 


393.26354 


Fe I 393.58213 


Ti II 


393.20172 


Fe I 


393.53133 





The ISS Ca K spectral band covers slightly more than 0.5 nm, centered on K3. 
Within this particular spectral range the ISS measures the profile of four neutral 
iron lines and a singly ionized titanium line. Table [2] lists these five spectral lines 
and their corresponding wavelength, taken from |Pierce and Breckinridge] (|1974p . 
These lines are used to determine the linear dispersion for the Ca II K range, as 
described below. 

4.1. Linear dispersion 

Each ISS measurement includes 1024 spectral samples covering the CaK spectral 
band. Although the ISS is an extremely stable instrument, mechanical hysteresis 
in the grating positioning mechanism may result in (sometimes significantly) 
changes in pixel position of spectra features. In order to address this issue, we 
first find the location of the K3 line core and then we use the relative distance 
of the other lines with respect to K3 to determine their positions. We use a 
quadratic fit and seven points around the line minima to determine the line core 
position of the five lines used for the calculation of the dispersion. 

Although the level 1 fits files provide the error in intensity for each pixel, we 
use unweighted polynomial fits throughout the data reduction process. The main 
reason for this choice is due to the large pixel-to-pixel variability of these errors, 
in some cases up to a factor 20 or more. When the uncertainty in intensity is 
included into the fit the fluctuations in the computed coefficients are significantly 
larger than in the case when the fit is unweighted, contrary to the very high sta- 
bility of the ISS. It also produces small pivot elements in the variance-covariance 
matrix associated with the observed points used in the fit of some spectra that 
results in a poor determinations of the parameters. 

We compute the formal 1-cr estimates of the returned parameters from the 
unweighted fit under the assumption that our model for the data is correct. 
This is certainly true for the high-quality data we are using. After the five line 
positions (in pixel units) are determined, a simple linear fit with the nominal 
wavelengths listed in Table [2] is used to calculate the linear dispersion. 

Figure [T] shows the comparison between the linear dispersion calculated using 
the new approach and the previous result from the ISS pipeline. A total of 
1678 observations has been used, that cover the period from December 2006 
to February 2011. Unlike the old reduction, all observations have now a well 
determined dispersion. The dispersion shows slight (~0.02%) annual variations 
that may be attributed to change in the index of refraction of air. A jump in 
the value of the dispersion occurring around June 2008 has been traced to a 
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Figure 1. Improved calculation of the linear dispersion for Ca II K observations taken during 
the period December, 2006 - February, 2011. The new determination of the five line core 
positions used to compute the dispersion clearly shows the stability of the algorithm over time 
(bottom plot). The error bars are not shown for clarity (see text). 



failure of a thermostat that controls the temperature in the spectrograph room. 
After the thermostat was replaced, the spectrograph had stabilized at a slightly 
different dispersion. However, one should notice that this change in dispersion 
is extremely small, about 0.1%. 

For purposes of clarity, we intentionally omitted the error bars in Figure [1] and 
all other figures showing the temporal behavior of the parameters computed from 
the Ca II K profiles. The errors on these parameters differ significantly depending 
on the strategy used to computed the parameters. With the approach described 
in this paper the errors are generally more uniform over time than the ones 
computed by the old ISS pipeline, in some cases smaller, and the corresponding 
error bars are typically shorter than the size of the points shown in the plots. 
We will discuss the errors in a separate section of this article. 

As part of this analysis we have investigated the temporal behavior of the core 
position and intensity for the five spectral lines used to determine the linear dis- 
persion. The core positions of each line shows no significant temporal trends but 
the fluctuations around the mean value, although very small, are quite different. 
These relative fluctuations vary from - 8 x 10~ 9 (Ti n 393.20172) to - 2 x 10~ 7 
(Fe I 393.58213). The core intensities of these spectral lines, on the other hand, 
show a tendency to increase over the time covered by this investigation. This 
increase is less than 0.005 in residual intensity, for all lines. 
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Figure 2. Comparison between the Ax3 and 1^3 time series computed from the old ISS 
pipeline (top) and the new and improved version (bottom). With the exception of very few 
data points, the impact of the new code in the computation of Ik3 is negligible. The error bars 
are not shown for clarity (see text). 



4.2. K3 wavelength position and intensity 

The wavelength position of K3, Xksi an d the corresponding core intensity Ik3 
can be determined using a quadratic fit over points centered around the pixel 
with minimum intensity. The choice of the number of pixels to be used in the fit 
plays a critical role in determining the quality of the result. The high spectral 
resolution of the ISS requires a relatively large number of pixels to properly 
characterize the core of K3, which is quite flat. For the new version of the 
reduction code we have used a total of 19 pixels. The result is a more stable 
determination of Xk 3 as shown in the bottom left panel of Figure [2] On the 
other hand, the impact on Ik3 is not significant. 

4.3. 1 A and 0.5 A emission index 

To calculate the 1 A and 0.5 A emission index we have used a new approach 
based on the Fourier components of the line profile. The central part of the Ca 
II K profile, between -0.6A to +0.6A from the center of K3, can be very well 
reconstructed using a limited number of Fourier components: 

NC 

I(x) = A + [ A i cos{2mx/N) + B t sin(2Trix/N)} , 
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Figure 3. Typical example of an observed ISS CaK profile (diamonds) reconstructed using 15 
Fourier components (solid line). The wavelength scale is measured with respect to the center 
of K3. 

where x = —N/2 + j (j = 0,N — 1), N is the number of spectral samples, 
and NC = 15 is the number of components used in reconstructing the profile. 
The coefficients Ai or Bi can be easily determined by multiplying both sides of 
the above equation by cos{2irkx / N) or sin(27rix / N) and summing over x, then 
making use of the property that sines and cosines are orthogonal. An example 
of an observed and reconstructed profile is shown in Figure [3J 

The major advantage of using this representation is that the above equation 
can be integrated analytically to compute the area under the line profile without 
the limitations in accuracy imposed by standard numerical integration schemes. 
Therefore, the Ca K EM emission parameter is given by: 

NC 

sin[7rz(a;2 — x-y)/N\ cos[7ri(x 2 + Xi)/N]}+ 

»=i 

JVC 

-2^{( J B l /i)sin[7ri(x 2 +x 1 )/N] cos[iri(x 2 - x^/N)}, 

i=l 

where x\ and x 2 are the integration limits. For example, in case of the 1 
A emission, x\ and x 2 are the values corresponding to -0.5 A and +0.5 A, 
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Figure 4. Comparison between CaK emission indices time series computed using a standard 
numerical integration scheme (trapezoid rule, top two plots) and the analytical formula de- 
scribed in this document (bottom two plots). The latest provides slightly higher values for the 
two emission indices. The error bars are not shown for clarity (see text). 



respectively. The resulting time series are shown in Figure 01 The new reduction 
indicates a slight increase in both the computed parameters. This is because 
the typical trapezoid rule applied for integration underestimates the actual area 
under the profile. 

4.4. K2 and Kl differences 

The wavelength positions of KIR (Xkir), K1V (Xkiv), K2R (\ K 2r), and K2V 
(^K2v) can be calculated using an approach similar to the one used to compute 
Xk3- However, in order to improve the quality of the result we included an 
additional step just before performing the quadratic fit. We take a 21-point 
interval centered around the minimum/maximum of K1/K2, and we fit this 
portion of the spectrum with a 6th order polynomial. Instead of using points 
from the observed spectrum, we use points from this polynomial for our 7-point 
quadratic fit around the positions of KIR, K1V, K2R, and K2V. The result is 
a more stable determination of the four wavelengths, as clearly shown in Figure 

El 

Once these four wavelengths have been determined, the calculation of the four 
corresponding intensities is performed as described in subsection 4.2. 
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Figure 5. Time series of wavelength separation of R and V emission K2 (left) and Kl (right). 
Shown is a comparison between the results obtained from the old ISS pipeline (top row) and 
the new calculation described in this article (bottom row). The error bars are not shown for 
clarity (see text). 



4.5. K2-K3 asymmetry and K2V/K3 ratio 

These two parameters are derived from K3, K2R and K2V intensities: 



K2V/K3 =I K 2v/Ik3 



K2-K3 asymmetry = (I K 2V - Ikz)/{Ik2R - ^3) 

These two parameter series marginally benefit from the new and improved 
version of the code, as shown in Figure [6l 



4.6. Wilson - Bappu index 
This index, WB, is defined as: 



WB = Iog[76.28(A a -Av)], 

where \r and Xy are the wavelength positions corresponding to (Ikir + 
Ik2r)/2 and (Ikiv + -Tkt2v)/2, respectively. In the original code these two 
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Figure 6. Ratio of K2V to K3 intensities and line asymmetry ratio of V and R parameter 
time series derived from the new ISS pipeline. The error bars are not shown for clarity (see 
text). 

positions were determined only to the rounded pixel. In this new version we 
use a proper interpolation to accurately estimate these two wavelengths. A 
comparison between these two versions of the code is shown in Figure [7] Note 
that the pixel-rounding problem clearly visible in the old time series has now 
been eliminated. 



5. Error estimate 

The original formulas used to calculate the errors associated with each parameter 
have not been changed in this new version of the pipeline. However, the error 
determination has greatly benefited from the improved estimation of the pa- 
rameters. Figure |5] illustrates a typical comparison between the formal 1-a error 
time series computed for Xk3 using the old ISS pipeline and the new improved 
calculation. The most important feature in this comparison is the much reduced 
scatter in the error determination from one observation to another. The much 
more homogeneous behavior in the error time series is achieved also for the other 
eight parameters. This result indicates that the error for all parameters has been 
been reduced significantly for many observed spectra as compared to the previous 
data reduction. This is particularly true for the measured wavelengths. Finally, 
we should mention that errors are not reported for the Wilson - Bappu index. 
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Figure 7. The calculated Wilson - Bappu index time series using the old ISS pipeline (top) 
and the new improved version (bottom). Errors are not computed for this parameter. 



6. Conclusions 

Time series plots shown in the previous section indicate significant improvements 
in parameters derived from Ca 11 K line profiles. For eight parameters, error bars 
are much smaller for many data points as compared with previous data reduction 
(no error bars are computed for Wilson-Bappu parameter). Most importantly, 
we see very little variation in the amplitude of error bars throughout the five 
year observing period. Together with nearly constant linear dispersion, this lack 
of variation in the amplitude of errors indicates the extreme stability of the 
Integrated Sunlight Spectrometer. The novel method of computing the emission 
indices based on the Fourier decomposition of the spectral line profile allows 
a more precise determination of these parameters. Comparison with widely 
adopted numerical integration schemes shows that these numerical integration 
methods may underestimate the emission measurements. 

Routine observations with the ISS were started at the end of 2006 at the 
declining phase of Solar Cycle 23. Although the scientific study of this data set 
is outside the scope of the current article, one can note clear cycle variations 
in several parameters of this spectral line with a minimum months earlier than 
the official sunspot minimum of December 2008. The K3 line core intensity, the 
two EM emission parameters, and the wavelength separation of V & R emission 
minimum show a minimum around August-September 2008 (± 40 days). The 
Wilson - Bappu index time series has an even earlier minimum, around the end of 
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Figure 8. Comparison between the formal 1-cr error time series for Aj^3 calculated using the 
old ISS pipeline (top) and the new improved version (bottom). 



Table 3. Cross-correlation matrix generated from the nine parameter time series. 
The indices correspond to: l=Line intensity (K3), 2=l-Angstrom emission index, 
3=Half- Angstrom emission index, 4=Ratio of K2V to K3 intensities, 5= Wavelength 
separation of V & R emission maximum, 6=Wavelength separation of V & R emission 
minimum, 7= Line asymmetry, 8= Wavelength for K3, and 9=Wilson-Bappu parameter. 





1 


2 


3 


4 


5 


6 


7 


8 


9 


1 


1.000 


0.992 


0.983 


-0.931 


-0.337 


0.655 


-0.548 


-0.629 


0.520 


2 


0.992 


1.000 


0.992 


-0.888 


-0.324 


0.693 


-0.582 


-0.654 


0.568 


3 


0.983 


0.990 


1.000 


-0.854 


-0.379 


0.761 


-0.638 


-0.709 


0.611 


4 


-0.931 


-0.888 


-0.854 


1.000 


0.249 


-0.397 


0.368 


0.433 


-0.295 


5 


-0.338 


-0.324 


-0.379 


0.250 


1.000 


-0.362 


0.214 


0.362 


-0.161 


6 


0.655 


0.693 


0.761 


-0.397 


-0.362 


1.000 


-0.820 


-0.809 


0.768 


7 


-0.548 


-0.582 


-0.638 


0.368 


0.214 


-0.820 


1.000 


0.795 


-0.694 


8 


-0.629 


-0.654 


-0.709 


0.433 


0.362 


-0.809 


0.795 


1.000 


-0.647 


9 


0.520 


0.568 


0.611 


-0.295 


-0.161 


0.768 


-0.694 


-0.647 


1.000 
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June 2008 (± 40 days). Such mismatchs between solar cycle minima in different 
indices has been reported in earlier papers (e.g., IKanel 12002)) . However, IKanel 
(2002| had found that Ca K-line intensity is in good agreement with sunspot 
number in timing of solar minimum, while our data show Ca II K3 minimum 
intensity preceding sunspot minimum. 

Finally, to better understand the relationship between pairs of parameters we 
have generated the cross-correlation matrix from the nine time series. This cross- 
correlation matrix is shown in Table[3] It appears that K3 intensity, 0.5-A and 1- 
A emission indices are highly correlated with each other. Wavelength separation 
of V & R emission maximum does not correlate well with any other parameters. 
K3 wavelength shows moderate correlation with most other parameters, and in 
fact, it varies inversely with solar cycle. 
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